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MODEUNG AT MORE THAN ONE LEVEL OF RESOLUTION 
Background of the Invention 
5 This application relates to representing computer models and more particularly 

to representing a geometry model at more than one level of resolution'. 

The accurate representation of subsurface topology can have a profound effect 
on the interpretation of a geoscience model. This is mainly due to the presence of material 
properties such as, for example, oil. For a more detailed introduction on the importance of 
10 topology, see United States Patent Application Serial No. 08/772,082, entitled MODELING 
GEOLOGICAL STRUCTURES AND PROPERTIES . 

Imagine, for example, the situation of two compartments 10a, 10b separated by 
a sealed fault 12 (a sealed fault is an impermeable membrane that does not permit fluids to 
pass), as shown in Fig. 1. Due to the sealed fault there can be no flow of fluid from 
15 compartment 10a to compartment 10b. If both compartments were to contain oil it would 
be necessary.to drill into both compartments to recover the oil. 

Now suppose the sealed fault 12 is punctured, as shown in Fig. 2. There is now 
a free flow of fluid between compartment 10a and compartment 10b so it would be possible 
to drill just into one of the compartments to extract the oil. 
2 0 Thus, having the correct topology in the geoscience model can have a profound 

effect on the finances of an oil-field development. 

Another important concept in interpreting geoscience models presented 
graphically on a computer screen is the concept of multiresolution analysis, whereby an 
analyst can view an area of interest in the geoscience model at different resolution levels. 
2 5 There are many techniques that have been developed for multiresolution analysis of 
surfaces. The work falls into two principal categories. The first is to use wavelets. The 
second is to use edge contraction and edge flipping, which is sometimes, called "topological 
editing". 

Wavelets have found wide acceptance in image processing and recently have 
30 found application in surface representations. Typically, wavelets are used in a subdivision 
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fashion. A typical subdivision scheme uses quaternary subdivision. For example, as shown 
in Fig. 3, a triangle 14 may be subdivided into four triangles 16a-d. In the example shown 
in Fig. 3» the retessellation conforms to the subdivision scheme. That is, the retessellation 
of the triangle into four triangles conforms with the quaternary subdivision scheme. It can 
5 be imagined, however, that if a different local retessellation of the triangle is performed, it 
may not be clear how to rebuild the subdivision scheme, since the refinement may not 
conform to the subdivision scheme. For example, if the triangle 14 is retessellated as in Fig. 
3b into triangles 17a and 17b, the retessellation does not conform to the quaternary 
subdivision scheme. There has been no work on the integration of wavelets and boundary 

10 representations. 

Topological editing, or editing a mesh using the topological operations of vertex 
removal, edge contraction and edge flipping, can be used to build multiresolution surfaces. 
Some mesh building techniques build a history of topological operations which permits 
progressive and partial loading of the surface, but it is not clear how this history is modified 

15 if a triangle is refined. There has been no work on the integration of topological editing and 
boundary representations. 

Summary of the Invention 
In general, in one aspect, the invention features a method, computer system and 
computer program for representing a first surface at multiple levels of resolution. The 

2 0 first surface comprises zero or more zero-cells, zero or more one-cells and one or more 
two cells. The method is implemented in a programmed computer comprising a 
processor, a memory, a persistent storage system, at least one input device, and at least 
one output device. The method and a model are stored on a computer-readable media and 
the method represents the model on one of the output devices. The method comprises 

2 5 partitioning the first surface with one or more boundaries, each level of resolution having 
a subset of the boundaries. 
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Implementations of the invention may include one or more of the following. The 
first surface may be partitioned into n\ nodes at resolution level-i using the level-i subset 
of boundaries. Each level-i+1 node may be associated with a unique level-i node. Each 
level-i node may be associated with the level-i+1 nodes associated to the node. Each 
5 level-i node may be associated with a subset of vertices that are critical at resolution level 
i. Assuming level d is the deepest level of resolution and the first surface is divided into 
simplices, each node at resolution level d may be designated a leaf node, each simplex 
may be associated to a unique leaf node, and each leaf node may have associated with it 
the simplices associated to that leaf node. 

10 A level-i node may have associated with it the list of simplices which is the union 

of all simplices associated with the level-i+1 nodes associated to the level-i node. The 
subset of boundaries for each node may be assigned to be the boundary of the union of the 
simplices associated with that node. The nodes may form an original tree and each node 
may be assigned a unique key. Each vertex in a leaf node may be assigned the key 

15 corresponding to that leaf node. 

The representation of a second surface may be stored in the computer-readable 
media, the second surface having nodes, leaf nodes, vertices, critical vertices and 
simplices. It may be determined which leaf nodes of the first surface intersect the leaf 
nodes of the second surface and the intersecting simplices from the first and second 

2 0 surfaces from the simplices associated to the intersecting leaf nodes. 



• 
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Each node except the leaf nodes may have a subtree, and the original tree may be 
split into new trees and each new tree may be associated with a new cell. The subtrees of 
the original tree which have no intersecting leaf nodes may be identified with one of the 
new cells. The simplices of the first surface may be split along the intersection curve. 
New simplices may be formed by tessellating the split simplices to respect the macro- 
topology of one-cells and zero-cells passing through the original simplices. A new tree 
may be built for each new cell and each new simplex may be assigned to the leaf node of 
the tree created for the new cell to which the new simplex belongs. For each leaf node of 
each new tree, each simplex in the original tree which is connected to a new simplex in 
the new tree leaf node and which lies in the same tree leaf node as the new simplex may 
be migrated. The neighbors of a tree node may be determined by finding all the keys of 
all the critical vertices in the node. The coarsest level node which is an ancestor of a key 
from the critical vertices in the migrated tree nodes and has not been split or migrated 
may be determined and that node may be migrated to the new tree. 

A complete node front of the tree of the first surface and a collection of vertices 
on a boundary of the first surface may be defined. A list of critical vertices from the tree 
nodes of the complete node front may be built. Those vertices identified to one- or zero- 
cell vertices may be removed from the list and all zero-cell vertices from the model which 
lie in the first surface and the defined collection of one-cell vertices may be added to the 
list. The collection of one-cell edges may be recorded. The surface may be tessellated to 
respect the list of vertices and the recorded one-cell edges. 

The subset of vertices on the boundary of the first surface which are also on the 
boundary of the second surface may be required to be the same as the subset of vertices 
on the boundary of the second surface which are also on the boundary of the first surface. 

A geometrical representation of the first surface may be maintained in a persistent 
storage. 
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A bounding box for each node may be stored on a persistent storage device, and 
for each node a list of critical vertices associated with that node may be stored on the 
persistent storage device. Storing the list of ' critical vertices may comprise storing a 
vertex descriptor for each vertex, storing a parameter value for each vertex, and storing an 
5 image value for each vertex. That portion of the first surface required may be loaded on 
demand from persistent storage and that portion of the first surface not required may be 
removed from memory. A tree node may be loaded from persistent storage on demand 
and removed from memory when it is no longer required. Simplices associated with a 
tree leaf node may be loaded from persistent storage on demand and removed from 
0 memory when no longer required. 

In general, in another aspect, the invention features a method, computer system 
and computer program for representing a first surface at multiple levels of resolution. 
The method is implemented in a programmed computer comprising a processor, a data 
storage system, at least one input device, and at least one output device. The method and 
5 a model are stored on a computer-readable media and the method represents the model on 
one of the output devices. The method comprises storing a grid representation of the first 
surface, the grid representation being made up of grid cells, forming a mesh 
representation of a portion of the first surface by triangulating a subset of the grid cells, 
and inserting the first surface into the model. 
0 Brief Description of the Drawings 

Figs. 1, 2, Sa, Sb, da^ and 35-41 are representations of models. 
Figs. 3a and 3b are representations of prior art multiresolution analysis schemes. 
Fig. 4 is a block diagram. 

Fig. 6b is a representation of the geometry of a model. 
5 Fig. 6c is a topology graph. 

Figs. 7a-7d illustrate simplices. 
Figs. 8a, 8b and 8c illustrate a simplicial complex, a simplicial complex with a - 
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hole and an illegal simplicial complex, respectively. 

Fig. 9 illustrates a non-manifold triangulated surface. 

Figs. 10a- 10c illustrate classification of one surface into another. 

Figs. 1 la-1 Id illustrate the need for coherent models to perform editing. 
5 Figs. 12a- 12c illustrate making a model coherent. 

Figs. 13a- 13c illustrate grid representations. 

Fig, 14 illustrates a mesh representation of a surface. 

Fig. 15 is a table of rendering times and rates. 

Figs. 16a-c illustrate cracking caused by approximation. 
1 0 Figs. 17a-b illustrate cracking caused by approximation. 

Figs. 18a-b illustrate bubbling caused by approximation. 

Figs. 19a-b illustrate bubbling caused by approximation. 

Figs. 20a-b illustrate the two possible triangulations of a grid cell. 

Fig. 21 is a graphical representation of a quadtree. 
1 5 Fig. 22 is a geometrical representation of a quadtree. 

Fig. 23 is a table showing the size of the different components of a quadtree. 

Figs. 24a-b illustrate the effect of removing a non-critical vertex through 
approximation. 

Figs. 25a-b illustrate the effect of removing a critical vertex through 
2 0 approximation. 

Fig. 26 illustrates intemal critical vertices. 
Fig. 27 illustrates external critical vertices. 

Figs. 28a-b, 29a-b illustrate approximation of a surface while respecting critical 

vertices. 

2 5 Figs. 30a-e and 33 are flow diagrams. 

Fig. 31 represents a quadtree with leaf nodes at a fixed depth. 

Fig. 32 illustrates parent critical, edge critical and sub-critical vertices. 
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Fig, 34 is a table listing common surfaces and their Euler Characteristics. 

Fig. 35 illustrates embedded disc obstructions that may be found in an external 

boundary. 

Fig. 36 iUustrates a method for removing a particular type of obstruction. 
5 Fig. 37 illustrates a method for determining the contour formed from the set of 

all vertices connected to a recently defined contour. 

Fig. 38 illustrates projecting an upper contour of a truncated square-based 
pyramid onto the interior of a square defined by the lower contour of the pyramid. 

Fig. 39 illustrates square based pyramids with contours meeting in a non- 
1 0 manifold manner. 

Fig. 40 illustrates contours being the boundary of an embedded disk. 
Fig. 41 illustrates pyramids formed with triangle caps and square bases. 

Description of the Preferred Embodiments 
15 . A geoscience environment 18, called GeoFrame, comprises an application 

framework 20 which comprises a shareable database 22, a geometry modeler (SGM) 24, a 
geometry query interface (GQI) 26, a solid modeler 28, an interactive geometric modeling 
library (IGM) 30 and a renderer (Visualize) 32, as shown in Fig. 4 and explained in United 
States Patent Application Serial No. 08/772,082, entided MODELING GEOLOGICAL 
2 0 STRUCTURES AND PROPERTffiS, incorporated by reference. 

An application accesses shareable database instance data using the Application 
Data Interface (ADI) 34. The ADI provides a programmatic interface to the shareable 
database 22. The geometry modeler (SGM) 24 describes the part of the shareable 
database 22 specific to 3D geometrical modeling- An application accesses SGM instance 
2 5 data using the GQI 26, 

The GQI provides a procedural interface to access the SGM. The GQI also 
provides a stable interface to a low-level non-geologic solid modeler 28 and adds extra 
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organization to the low-level geometric representations to aid geological understanding. 

The GQI is implemented on top of a low-level geometric solid modeler 28, such 
as SHAPESTM from XOX. 

The Interactive Geometric Modeler (IGM) 30 provides visualization services on 
5 behalf of and a graphical interface to the GQI. An application normally receives human 
instruction through the IGM. 

The IGM is implemented on top of a low-level rendering engine 28, such as a 
combination of OPENINVENTOR and OPENGL. 

The Common Model Builder 20 is a toolkit that brings together the 3D 
10 geometric modeling components. It integrates data management services, geometric 
modeling and graphical interfaces for human input. The Common Model Builder provides 
an abstract framework upon which a 3D geoscience application can be built. 

The GQI clearly distinguishes between the notions of topology and geometry. 
Broadly, topology refers to the connectivity between components in the model and generally 
15 refers to "macro-topology", described below. Geometry refers to the actual point-set 
representation of a particular component, for example. 

Macro-topology carries the relationships between the major topological 
components of a model. The major topological components refer to points, edges (curves), 
faces (surfaces), and subvolumes. For example, macro-topology would answer a question 
20 such as "which surfaces bound a particular volume?" or "which surfaces lie within a 
particular volume?". Macro-topology does not need to consider the underlying geometry of 
surfaces and can work with many different types of geometrical representations, for 
example NURBs and meshes. Macro-topology is implemented in the SHAPES solid 
modeler 28. 

25 In the terminology of the GQI, macro-topology is represented by cells. A cell is 

a path-connected subset of Euclidean space of a fixed dimension. Path-connected means 
any two points in the cell can be connected by a path in the cell. The dimension refers to 
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the dimension of the geometry of the cell, a 0-cell is a point, a 1-cell is an edge (curve), a 2- 
cell is a face (surface) and a 3-cell is a volume. For example, areas 36 and 38 are distinct 
cells separated by fault 40, as shown in Fig. 5a. Curve 42 illustrates the path-connected 
character of cell 38. Similarly, as shown in Fig. 5b, the line segments 44 between 
5 intersection points and the intersection points 46 themselves are cells. 

The topology of a geometric model is the set of all cell-cell connectivity 
relationships. The cell-cell connectivity relationships can be represented in a graph in 
which the arcs of the graphs represent connectivity and the nodes represent cells, as shown 
in Figs. 6a-c. In a cellular model, cells of dimension n are connected to boundary cells of 
10 dimension n-1 and vice versa. For example, in the topology of a box 48, shown in Fig. 6a, 
the area cell 50 is connected to its four bounding edge cells 52a-d. A single cell, can act 
both as a boundary and as a region. For example, a surface can bound a subvolume, but can 
itself be bounded by a set of curves. The 1-cells 52a-d, in Fig. 6b, are both subregions 
(bounded by zero-cells 54a-d) and boundaries (of area cell 50). These relationships can be 
15 represented graphically, as shown m Fig. 6c. Node 56, corresponding to area cell 56. is 
connected to nodes 58a-d, representing 1-cells 52a-d, respectively. The connection between 
the nodes is represented by the arcs between them. 1-cells 58a-d are connected to zero-cells 
60a and 60d, 60a and 60b, 60b and 60c, and 60c and 60d, respectively. Cells which are 
contained in higher-dimensional cells but do not split them (such as fault 40 in Fig. 5a) are 

2 0 said to be "embedded". 

Geometry is the point-set description of a topological object. For example, a 
geometrical query would be "which is the closest point on the surface to a given point?". 

All geometry to be considered here is piecewise linear. For convenience, a brief 
review of the definition and topological properties of a simplex, simplicial complex, and 
25 triangulated surface or mesh is included below. In short, the piecewise linear geometry of 
an 1-cell or edge is a polyline, for a 2-cell or face it is a triangular mesh and for a 3-cell it is 
a tetrahedral mesh. 
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Let ao,...,ar be a collection of r+1 linearly independent points in R" with r<n. 
The r-dimensional simplex a = (ao,...,ar) is the set of all points in R" such that 

<T = L € /? " I X = X . a, , > 0, X r, = 1 ► 

I 1=0 1=0 

The r-tuple (to,...,tr) is the barycentric coordinate of x. The [a\] are the simplex's comer 
5 points. 

A 0-simplex is a point 62, illustrated in Fig. 7a, and a 1-simplex is the line 
segment 64 joining, for example, ao to ai, as shown in Fig. 7b. A 2-simplex is a filled 
triangle in 3-space 66, the triangle being defined by the three comer points, as shown in Fig. 
7c. Similarly, a 3-simplex is a filled tetrahedron 68 defined by the four comer points, as 
1 0 shown in Fig. 7d. 

Let a = (ao,. .,ar) be an r-simplex and let q be an integer with 0 < q < r. Now 
choose q+1 distinct points from ao,....ar, say aio,...,aiq, these q points define a q-simplex 
which is called a q-face of a. 

An n-simplex has n+1 n-l-faces. These faces are exactly the faces expected, that 
15 is the 0-faces of a 1-simplex are the end points, the 1 -faces of a 2-simplex are the bounding 
lines, the 2-faces of a 3-simplex are the triangular faces of the tetrahedron. 

A simplicial complex, K, is a finite collection of simplices in some R" satisfying 

1 . If a e K, then all the faces of a belong to K 

2. If a, T G K then either a n x = null or a o i: is a common face of c 
20 andx. 

Examples of legal simplicial complexes include a simplicial complex with 
vertices, e.g. 70, and simplices, e.g. 72, as shown in Fig. 8a, and a simplicial complex with a 
hole 74, as shown in Fig. 8b. Overlapping simplices, e.g. 76, create an illegal simplicial 
complex, as shown in Fig. 8c 
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The "dimension" of a simplicial complex is -1 if K = null, and the maximum 
dimension of the simplices of K otherwise. 

Given a collection of simplices C from the simplicial complex K, the collection 
C defines another simplicial complex by considering all the faces of C: 
SimpCompiC) = {(J 6 K\aisa faceof asimplexinC^ 



Let K be a simplicial complex and define: 



risan s — simplex 

T is the face of exactly r distinct 5 + 1 simplices 
T is not the face of any 5 + 2 simplices 



and 



SI iK) = SimpCompiA'^ (K)) 

10 A pseudomanifold, M, is a simplicial complex such that 

1. M is homogeneoiisly n-dimensional. That is, every simplex of M is a 
face of a n-simpiex of M. 

2. Every (n- 1 )-simplex of M is a face of at most two n-simplices. 

3. If a and tf are two distinct n-simplices of M, then there exists a 
15 sequence ai,....ak of n-simplices in M, such that ai = a, Ok = a' and oi meets 

Qi+i in a (n-l)-face for 1 < i ^ k-1. 

A triangulated surface or mesh is a two-dimensional pseudo-manifold. The 
presence of a singular vertex 78, as shown in Fig. 9, means that the triangulated surface 80 
is not a manifold, but it is a pseudo-manifold. 
2 0 For an n-homogenous simplicial complex K, i.e. every simplex is a face of some 

n-simplex, the following can be seen to hold: 

S^iK)^S^(K) = .,.=^Sl,(K) = null 
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5 




an homogenous simplicial complex of dimension n-1. However, the boundary of a . pseudo- 
manifold is not necessarily a pseudo-manifold. 

Denote by Simps(v) the set of 2-simplices which connect to the vertex v in a 
triangle mesh. 



geometry engine. The SHAPES geometry engine uses micro-topology for this purpose. 



representation is an integral part of the surface representation. It is sufficient for the mesh 
15 representation to support basic navigation queries, such as, "which simplices use a vertex?" 
and "which vertices are used by a simplex?". Using these two queries it is possible to 
navigate from a simplex to a vertex and then from this vertex back to another simplex. In 
this way it is possible to determine connected components in a model. Micro-topology 
supports these types of queries. 
2 0 An additional requirement, which is imposed by the SHAPES geometry engine, 

is that all intersection curves must be supported by micro-topology. 



vertices lie on the boundary of a surface. This information is maintained by the SHAPES 
geometry engine and permits efficient queries identifying whether a vertex is identified to a 
2 5 1-cell vertex or a 0-cell vertex. 



10 



Denote by Verts(S) the set of vertices forming the comers of the 2-simplex S. 
Any surface representation must be able to pass geometry to the underlying 



Micro-topology is an example of a mesh representation. The mesh 



For the purposes of the surface representation it is necessary to know which 
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An additional requirement imposed by the SHAPES geometry engine is that all 
vertices which are identified with a 1-cell or 0-cell vertex and all boundary simplices are 
retained in core memory at aU times. 

One of the algorithms provided by the GQI is classification. Given two cells, A 
and B, classification subdivides the respective point sets into an inside part, an outside part, 
and a part on the boundary of the other. In set theoretic terms, the union of A and B breaks 
into three disjoint components, the part of A which is not in B (ANB), the part of B which is 
not in A (B\A) and the intersection of A and B: 

A o B = A\B B\A u ( A o B ) 
A\BoB\A = null set 
AVB n (A o B) = null set 
B\A n (A o B) = null set 
When performing the classification, the GQI identifies the connected components or cells in 
the model, connects these cells along their shared boundaries, throws away the original cell 
definitions and builds new cell definitions for the connected components. The outcome is 
an Irregular Space Partition (ISP), 

For example, as shown in Figure 10a, if an earth model 82 initially comprising a 
volume 84 and two horizons 86 and 88 has inserted into it a fault 90, shown in Fig. 10b, the 
classification proceeds as follows. 

Compute the intersections of shapes and generate cells representing the 
intersection geometries. In Figure 10c, the intersection geometry is represented by the 
heavy dot intersection points 92. 

Split cells are sub-divided by lower-dimensional cells. In Fig. 10c the horizons 
86 and 88 intersect fault 90. The fault 90 ceases to exist and is replaced by cells 94, 96 and 
98. Similarly, the combination of horizons 86 and 88 and cell 96 splits area 84 in two. 
Thus, area 84 ceases to exist and is replaced by areas 100 and 102. 
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Classification is a macro-topological operation which is implemented using 
micro-topological functionality. All surfaces that provide micro-topology must be able to 
support the micro-topological interface required by the macro-topological interface to 
perform classification. 

5 "Coherency" occurs if the geometry of the cells agrees in all dimensions. The 

need for coherency is illustrated, for example, in Fig. 1 la-d, where a volume of interest 104 
is classified with a saltdome 106 and classified with a horizon 108a-c. That part of the 
horizon 108b which is interior to the saltdome and geologically would not be present is then 
removed. 

10 To remove the interior surface 108b requires breaking up the horizon into two 

components, an exterior part 108a, 108c and an interior part 108b. To do this requires two 
distinct geometrical parts, the interior part and the exterior part. Once the geometry has 
been identified it is then a matter of removing the geometry and recomputing the macro- 
topology. This operation can be performed if the model is coherent. 

15 An ISP is coherent if the geometry of the cells agrees at all dimensions (i.e., if it 

is a homogenous simplicial complex, as defined above). Effectively, this means the 0-cell 
0-simplices and 1-cell 1-simplices are faces of the 2-cell 2-simpiices. 

Retuming to Figs, lla-d, to enable the removal of the interior surface 108b 
requires a geometrical representation of the surfaces which is coherent. In Figs. 12a-c, a 

2 0 cross-sectional view of an example of classification followed by coherency is given. A cell 
110 consists of two simplices 112, 114, as shown in Fig. 12a. Cells 112 and 114 are 
coherent because their geometry agrees at all dimensions. Volume 110 is coherent. Cell 
1 10 is classified against surface 116 which results in two new cells, cell a = {118, 122} and 
cell b = { 120, 124 } . Cells a and b are not coherent because their boundaries are not faces of 

2 5 simplices. 

To correct this, simplices 118, 120, 122, and 124 are split to respect the 
boundaries of cells a and b, into simplices 126a-d, 128a-c, 130a-c and 132a-d, respectively. 
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This makes ceDs a and d coherent. 

There are two stages to making a model coherent. The first requires a re- 
tessellation of the surface in the neighborhood of the intersection curves in order to ensure 
the model is coherent. The second is to perform a migration, which involves collecting 
5 together the simplices into connected components or cells and identifying each connected 
component with an appropriate cell in the ISP. 

The retessellation can introduce tolerancing issues, because in performing the 
retessellation it is possible that degenerate triangles are produced. For example, a 
degenerate triangle may have a poor aspect ratio, or perhaps a small area. It is the 
1 0 responsibility of the underlying geometry engine to handle the degenerate triangles, which is 
typically done by merging degenerate triangles into neighboring non-degenerate triangles. 

In the GQI, it is necessary for the surface representation to be "parameterized". 
This means each surface has an "image space" and a "domain space". The image space has 
to be the same for all surfaces and is typically three-dimensional Euclidean space. The 
1 5 domain space for a surface is typically a rectangle in the plane and as classification proceeds 
the domain is restricted to subsets of the rectangle. 

There are two discrete surface representations which are of interest, grid and mesh. 
In any surface representation there are two important aspects, the topology and the 
geometry. Topology is used in its generic sense and refers to how the surface is connected, 
20 for example, which triangles connect to which triangles. The geometry specifies the actual 
position of the surface. The primary distinction between grid and mesh is grid represents 
topology implicitly, mesh represents topology explicitly. 



topology is represented implicitly. They differ in how the geometry is represented. The 
25 topology is represented implicitly as two integer extents giving the number of grid cells in 
each direction. 



There are several types of grid. They share a common feature which is that 
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The most compact fomi of grid is a "regular grid", as illustrated in Fig. 13a. For 
a regular grid it is only necessary to store an origin, step values for the grid points, the 
number of grid points and the height values for each grid point. A regular grid has a 
number of drawbacks. In particular, all of the grid cells are of fixed size and only height 
5 fields can be represented. 

A slightly more general form of grid is a "rectilinear grid". For a rectilinear grid 
the grid cell sizes can vary along each axis, as illustrated in Fig. 13b. As with regular grids, 
however, only height fields can be represented. 

The most general type of grid is a "structured grid", as illustrated in Fig. 13c. As 
10 with regular and recdiinear grids, the topology is represented implicitly by two integer 
extents giving the number of grid cells in each direction. The geometry is represented 
explicitly by maintaining a three dimensional point coordinate for each grid point. A 
structured grid differs from regular and recdiinear grids because it is possible to represent 
multi-valued height fields. It is also more adaptive than regular and rectilinear grids 
5 allowing the grid cell size to vary across the whole grid. 

Structured grids have the following characteristics: (1) they are compact; 
(2) they can represent multi-valued height fields; (3) topology is represented implicidy in 
that grid cell neighbors are given by increments and decrements of an indexing function; 
(4) the grid index is a natural parameterization; and (5) they are difficult to edit. 
0 The major drawback of all grid representations is the inability to topologically 

edit the grid. For example, it is easy to move a vertex in a structured grid by replacing the 
coordinates. It is difficult to insert a new vertex into the grid, which would require 
regenerating the grid indices for the surface. 

The greatest advantage of a mesh representation is the ability to represent 
5 irregular geometries. A mesh 142 with an irregular hole 144 can be represented very 
simply, as shown in Fig. 14. Where a surface is rapidly changing the mesh can be very fine, 
and in large flat regions the mesh can be very coarse. It is also very easy to edit a mesh. For 
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example, if a new vertex needs to be inserted, it is possible to retriangulate the surface in the 
neighborhood of the vertex. 

Mesh has the following characteristics: (1) triangles are of variable size; (2) they 
can represent multi-valued height fields; (3) topology is represented explicitly; (4) locally 
5 editable, for exaniple refinement and coarsening; (5) not necessarily parameterized. 

Meshes can represent more general surfaces than grids. Compared to a grid, 
however, a mesh incurs a memory and performance cost because the topology has to be 
represented explicitly. However, because a mesh is irregular it can easily model multi- 
valued surfaces and surfaces which are rapidly changing in one area and flat in another. 
10 When creating any surface representation, the sampling size can have a 

profound effect on memory usage. For example, if the samphng along each axis of a grid is 
doubled, the memory usage will be quadraplcd. This suggests the need to selectively load 
portions of the model in core memory. 

Furthermore, it is easy to build a geometric model which can reside in core 
15 memory, but which overwhelms even the most powerful graphics hardware. This suggests 
the need to be able to decimate or subsample the surfaces in an efficient manner for 
rendering. 

Consider the following scenario. A geologist wants to create and edit on a 
workstation a 3D geometric model made up of meshes containing on the order of 250,000 

2 0 simplices per surface and to have on the order of 50 surfaces. Fig. 15 records the rendering 
time in seconds on an SGI INDIGO^ EXTREME with a GUI-Extreme graphics card and 
128 Mbytes core memory for a set of N x N rectangular grids for a range of N. Each N x N 
grid contains 2(N-1)^ simplices. A grid*s nodes are evenly spaced in two dimensions, e.g. x 
and y axes, but not in the third dimension, e.g. z. An optimization tool by Open Inventor 

2 5 (ivfix) was used to optimize the output. 

In this example, the rendering rate varies with grid size, but is roughly 200,000 
simplices per second. For good interactive response, rendering times on the order of 0.1 
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second are required, because a geologist-user wants to be able to u-avel interactively around 
the model. In the example of an SGI INDIGO^ EXTREME, this speed is achieved for a 160 
X 160 grid. In general, this figure will be highly dependent on the particular graphics 
hardware installed on a particular machine. 
5 Increasing the power of the graphics hardware can significantly improve the 

rendering performance, but even this is of limited utility. Imagine a window on a computer 
of 500x500 pixels and hardware that is capable of rendering 50,000 triangles in 0.1 seconds. 
The total number of pixels is 250,000 that gives an average area of 5 pixels for each 
triangle. If the power of the rendering hardware is increased five-fold, a triangle will have 

10 an average area of 1 pixel. Beyond this point, increasing the power of the graphics 
hardware, and hence permitting more triangles to be rendered, will not significantly 
improve the quality of the screen image. 

Real world surfaces are typically non-uniform. In one region a surface may be 
varying rapidly, in another it may be almost flat. This suggests the need for allowing the 

15 user to selectively choose different levels of detail in different regions of the model. 

A geologist uses a 3D geometrical representation much like a microscope, in the 
sense that she wants to zoom in and zoom out in relatively small regions in a controlled 
manner. This usage is similar to that supported by terrain visualizadon implementations. 

For visueilization it would beneficial to draw parts of the model in the far field of 

2 0 view in a coarse representation, the parts of the model in the near field of view in detail. 

After classification and coherency procedures have been performed, a valid ISP 
has been built. Typically, the model is too detailed to be viewed interactively on hardware 
that exists at present so the amount of detail from the model is reduced or "decimated". To 
avoid unnecessary eutifacts in the model it is important that the decimated model respect the 

2 5 topology of the underlying ISP, as discussed above. This means all features of the boundary 
representation must be represented. The two effects that must be avoided are cracking and 
bubbling as discussed below. 
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There are two basic changes that can happen in the topology of a model. These 
two changes can happen in two places in the topology, the first is within a surface itself, the 
second in the boundary representation. 

Cracking is often seen as a hole where there was no hole before, as shown in 
5 Figs. 16a-c. Two approximations are shown of a half-cylinder 146, shown in Fig. 16a. The 
first approximation 148, illustrated in Fig. 16b, exhibits a hole 150 in the surface, while the 
second approximation 152, shown in Fig. 16c, does not. 

Cracking can also occur in the boundary representation of the BP, as shown in 
Figs. 17a-b. An ISP 154, shown in Fig. 17a, has an intersection curve 156 splitting a 
10 surface into two cells 158 and 160. When the ISP is decimated as shown in Fig. 17b (with 
cell 158 being approximated, but cell 160 not being approximated), the intersection curve 
on cell 158 has been approximated by a straight line 162, the intersection curve on cell 160 
has not been approximated. Due to the different approximations, the boundaries of cells 
158 and 160 are no longer coincident and a hole 164 appears. 
1 5 Bubbling is where intersections occur where there were none before decimation. 

As with cracking this can occur within the surface and within the boundary representation. 
Bubbling within a surface can occur where a non-self-intersecting surface becomes self 
intersecting, as shown in Fig. 18a-b. In Fig. 18a, a cross-sectional view is given of a surface 
166. Vertex 168 is dropped which causes the surface to become self-intersecting at point 
20 170, as shown in Fig. 18b, introducing a bubble 172. 

Bubbling can occur in a similar fashion to the cracking that occurred in Figs. 
17a-b, as shown in Figs. 19a-b. Fig. 19a shows a cross-section of an ISP 174 comprising of 
two surfaces 176 and 178. In a decimated version of the ISP 180, shown in Fig. 19b, surface 
176 is decimated by dropping vertex 182, while surface 178 remains unchanged. This 
25 introduces new intersections 184 and 186 between surfaces 176 and 178 that did not 
intersect in the original ISP 174, producing bubble 188. 

A new hierarchical surface representation is able to represent a surface that is 
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initially imported as a structured grid and is then edited by the use of classification and 
coherency. The representation supports the micro-topological interface required by the 
SHAPES geometry engine for a surface to be used in classification. The surface also 
supports adaptive decimation algorithms that prevent cracking and bubbling. The 
5 architecture is a hybrid-grid mesh. This means wherever possible a grid is used, but, when 
necessary, mesh regions overlaying the grid are used. This has many advantages, for 
example, being able to tune algorithms to make use of a grid representation when the 
surface is a grid. But, when greater flexibility is required, the grid representation can be 
exchanged for a mesh representation. 

10 Typically, surfaces are imported as grids. However, the geometry engine can 

only perform geometric calculations with mesh. The hybrid grid-mesh representation neatly 
solves this problem, because areas of the grid can be dynamically converted to mesh thereby 
supporting general topological and geometrical editing. But, the efficient grid 
representation can be used where such general editing is not required. Effectively, the 

1 5 underlying geometry engine is fooled into believing the surface is a mesh. Furthermore, 
since it is possible to identify which areas are grid and which are mesh, algorithms can be 
optimized to make use of the grid structure whenever possible. An example of tuning when 
the surface is a grid is encoding the particular triangulation of the grid. A grid cell has two 
possible triangulations 190 and 192, as shown in Figs. 20a-b. 

2 0 The triangulation of the grid can be stored in a bit vector, each bit representing 

the chosen triangulation of a particular grid cell. In this way it is not necessary to maintain 
the triangulation of the grid explicitly as simplices. 

hi general, the grid can be thought of as being the background and the mesh the 
foreground. Since the surface maintains its grid representation in the background it is 

2 5 possible to discard the mesh foreground at any time and recover the grid background. At 
any time a portion of the grid can be turned to mesh by reading the triangulation bit-vector 
and dynamically building simplices. 
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The hybrid architecture maintains the flexibility to provide the irregular 
refinement of a grid. This is important, for irregular refinement is required for efficient 
classification and coherency algorithms. 

For the purposes of the implementation, a quadtree has been used to provide a 
5 multiresolution hierarchy. A quadtree was chosen because of its geometrical relationship to 
sub-sampling in grids. 

A quadtree is a tree with nodes, each of which has four children, except for the 
leaf nodes, which have no children. Any quadtree can be drawn graphically, as shown in 
Fig. 21, or geometrically, as shown in Fig. 22. Every node of the tree has a unique depth 
1 0 and can be assigned a unique key. 

The key of a quadtree node is chosen to provide the following fiinctionality. 

• A compact and efficient way to dereference quadtree nodes. 

• A linear ordering for the quadtree nodes. 

• An efficient way to compute the depth of the key. 

15 • An efficient way to compute whether a key is an ancestor of another key. 

• An efficient way to compute whether a key is a descendant of another 
key. 

• A way to compute the ancestor keys of a key. 

• A way to compute the descendant keys of a key. 

2 0 Efficient in this context means a small number (typically <5) of bitwise Booleans together 
with bit shifts and arithmetic operations. The implementation of the quadtree key can be 
done by using pairs of bits to hierarchically identify the child at each level and representing 
the depth as an integer. Let ceil(x) be the smallest integer greater than or equal to x. Using 
the implementation of the quadtree key described above, the number of bits required to 

2 5 represent a tree of depth d is given by 2d+ceil(iog2(d+l)), where the factor 2d is the 
mantissa (2 bits are required to identify each child at each depth), and the factor 
ceil(log2(d+l)) is the number of bits required to represent the integer value d. Thus a 32 bit 
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key can represent a tree of depth 14, and as shall be seen this is more than adequate for 
current needs. 

The mantissa is implemented using bit interleaving. 

In Fig. 23 the size of the different components of a quadtree are given for 
5 different depths of the quadtree. It can be seen a tree of depth 10 can support a grid of size 
1000 by 1000, furthermore, the key is comfortably maintained in 32 bits. 
The following definitions will be used: 

• The number of elements (cardinality) in a collection C will be denoted 
by Card(C). 

10 • The key of a quadtree node N will be denoted by Key(N). 

• The quadtree node of a key K will be denoted by Node(K). 

The depth of the key K will be denoted by Depth(K). Without loss of 
generality, the root key has depth 0. 

• The ancestor key at depth i of the key K will be denoted by Ancestori(K). 
1 5 The function is defined for i < Depth(K) with AncestorDepth(K)(K) = K. 

• Denote by Ancestorsi(C) the collection of ancestor keys of the collection 
of keys C: 

Ancestors. {C)- \J Ancestor.{K) 
Note there is a one-to-one correspondence between nodes and keys: 

Node{Key{N)) = Key(NodeiK)) = K 
For simplicity, and since the original surface is a structured grid, the quadtree 

2 0 leaf nodes are assumed to be all at the same fixed depth. It is possible for branches of the 

quadtree to be empty. To ensure unnecessary navigation across empty parts of the tree, if a 

node is present in the tree it must have a non-empty leaf node in its descendants. 

Conceptually, all surfaces are represented as mesh. However, it is not necessary 

for the mesh to be completely built, in fact, it is possible to dynamically build the mesh as 
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and when required provided sufficient topology is maintained to support geometrical and 
topological algorithms. Once an area of the surface is marked as mesh and simplices have 
been built for this region, more general, topological editing can be performed in this region, 
for example, refinement. 
5 Following this paradigm, simplices are conceptually assigned to quadtree leaf 

nodes in a regular manner. As discussed above, this can be easily done for a grid. The 
quadtree leaf node maintains a flag that signifies whether its simplices have been built or 
not. If asked for its list of simplices the leaf node can return the list of simplices if they 
have been built, or build them dynamically and return the list. In this way a simplex is 

1 0 assigned to a unique quadtree leaf node. 

The geometrical representation of a quadtree mirrors the regular structure of a 
grid and hence it is natural to assign each grid cell in a regular manner to a unique quadtree 
leaf node. For example, the geometric tiling of a quadtree described above can be used and 
the grid cell may be assigned to the quadtree leaf node that contains its lower left comer. 

15 When a grid cell is triangulated, it contains a pair of simplices. These simplices are 
assigned to the quadtree leaf node of their grid cell. 

The quadtree leaf node that contains the simplex S will be denoted by Leaf(S). 
Having made the assignment of simplices to quadtree leaf nodes, the quadtree 
hierarchy gives an equivalence relation for each depth i of the quadtree, by assigning 

2 0 simplices to their ancestor node at level i. 

If S and S' are simplices then 5 - 5' => Ancestorf{Leaf (S)) = Ancestor.(LeafiS')) 

In other words, each quadtree node is assigned the simplices of its descendants and we 
define for a quadtree node: 

SimpsiN)^{S\S G N} 
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As has been explained, each quadtree node has been assigned a collection of 
simplices. The quadtree node inherits a boundary firom its collection of simplices, i.e., the 
boundary of the simplices. This boundary is a list of vertices and edges connecting the 
vertices. Furthermore, at a fixed depth in the tree, the boundaries of all the quadtree nodes 
5 at the chosen depth form a graph. 

The graph of the boundaries of the quadtree node will be a particularly detailed 
object for it takes its edges from the simplices in the full-resolution surface. It would be 
preferable to approximate the boundary of the quadtree node with a reduced set of edges 
and hence a reduced set of vertices. For example, as shown in Fig. 24a, it can be seen that 
1 0 vertex 194 can be dropped and the edge 196 which it defines. The edge can be straightened, 
as shown in Fig. 24b, without changing the basic structure of the graph (more precisely, the 
topology of the graph has not changed). 

In contrast, if vertex 198, shown in Fig. 25a, is collapsed to vertex 200, the 
valence of vertex 200 changes from three to four and hence the topology of the graph is 
1 5 modified. It can be seen that the only vertices which can be dropped and still preserve the 
topology of the graph are the vertices of valence two (such as vertex 194 in Fig. 24a). This 
can be made precise by introducing the notions of homeomorphic graphs. Graphs are well- 
known topological objects in mathematics and many techniques have been developed to 
study them. 

20 From the description above, there is a collection of vertices that cannot be 

removed without changing the topology of the graph. These vertices are of interest as they 
encode the topology of the surface. However, it can be difficult to compute this collection 
of vertices. The following discussion describes a computationally cheap means to find a 
collection of vertices, called "critical vertices", which include the vertices described above. 

2 5 The critical vertices are the crucial component when describing the topology of 

the surface. They enable navigation to be performed within the tree without requiring the 
complete tree to be loaded in memory. Moreover, they avoid creating unnecessary mesh 
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regions to describe the geometry of the surface. The navigation is generally an important 
part of a surface description but is essential in the operation of making coherent 

The critical vertices are also crucial in the decimation stage because they are the 
vertices that will appear in the decimated model. 
5 There are two classes of critical vertices: "internal" critical vertices and 

"external" critical vertices. The internal critical vertices are present to provide topological 
connectivity in the interior of the surface. The extemal critical vertices provide topological 
connectivity around the boundary of the surface and along cracks in the surface. 

Again, the driving example is the sub-sampling of a regular grid. The critical 
1 0 vertices mirror sub-sampling in a regular grid. 

A vertex is an internal critical vertex if the vertex is in the interior of the surface 
and it can not be removed from the graph of edges without changing the topology of the 
graph. A vertex is critical at depth i if it is at the intersection of three or more quadtree 
nodes at depth i. For the regular partition of a grid, the critical vertices are the interior 
1 5 vertices of the sub-sampled grid, as can be seen in Fig. 26. 

The internal critical vertices do not include the boundary vertices. To include 
these requires the macro-topology of the 2-cell to be used and gives the notion of an 
extemal critical vertex. 

A vertex is an extemal critical vertex at depth i if it is: 
20 ♦A vertex which is identified with a 0-cell. 

• A vertex which is identified with a 1-cell vertex and lies at the boundary 
of two or more quadtree nodes at depth i. 
The extemal critical vertices permit the boundaries to be included, as can be seen in Fig. 27. 

The collection of critical vertices is the union of the internal critical vertices and 
2 5 the extemal critical vertices. One can see from the definitions if a vertex, v, is critical at 
depth d then the vertex is critical at all depths greater than d. 

The "depth" of a critical vertex is the depth at which the vertex first becomes 
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criticai. If a vertex is never critical in the tree, it will first appear in the vertices of the 
simplices of a leaf node and its depth is given by the depth of the tree plus 1 . The critical 
vertices represent an approximation of the surface and at greater depths of the quadtree the 
approximation improves. Figure 28a shows an original surface with a hole comprising two 
5 1 -cells 198 and 200 and a crack comprising one 1-cell 202 with a quadtree overlaid. At this 
depth, there is only internal critical vertex 204 (occurring at the intersection of three of more 
quadtree nodes). There are external critical vertices 206, 208, 210. 212, 214, 216, 218, and 
220 identified with a 0-cell. There are 9 external critical vertices 222, 224, 226, 228, 230, 
232, 234, 236 and 238 identified with a 1-cell vertex of the surfaces and lying at the 

1 0 boundary of two or more quadtree nodes. Fig. 28b shows the approximation of the surface, 
hole and crack at this level of decimation. Curves between critical vertices are replaced by 
straight lines. For example, the portion of curve 198 between critical vertex 206 and critical 
vertex 228 is replaced by line 239. 

Fig. 29a shows the same original surface overlaid with a quadtree at a depth one 

1 5 greater than that shown in Fig. 28a. The number of quadtree nodes quadruples from four to 
sixteen. Consequently, seven interior critical nodes 240, 242, 246, 248, 250, 252 and 254 
are added. Further, 1-cell external critical vertices at the boundary of two or more quadtree 
nodes 256, 258, 260, 262, 264, 266, 268, 270, 272, 274, 276, 278, 280, 282, 284 and 286 
are added. Fig. 29b shows the approximation of the surface, hole and crack at this level of 

2 0 decimation. Again, curves between critical vertices are replaced by straight lines. 

As can be seen, assigning critical vertices to a quadtree node allows a decimated 
view of a surface to be built. For the quadtree node, N, let 

CriticalVertsiN)=^{v\veVertsiS) for some S e Simps(N)anddepthiv) <depthiKeyiN))} 

In its most abstract form, a 2-simplex is formed from three vertices. In a mesh, 
25 the connectivity of the surface can be represented by sharing the same vertex across 
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different simplices and maintaining within a vertex references to the simplices which 
reference the vertex. So, given a simplex, one can navigate to one of its vertices and from 
the vertex can navigate back to another simplex. In this manner one can travel around the 
simplices and vertices in a surface. 
5 In the hierarchical surface representation, connectivity is built in a different 

manner. Within each critical vertex, instead of storing the simplices which reference the 
vertex, the keys of the quadtree leaf nodes which have the vertex as a critical vertex are 
stored. This list of keys is called the vertex descriptor. In this manner, it is not necessary to 
instantiate simplices to be able to navigate around the surface. Navigation can be performed 

10 from a quadtree node by using the quadtree leaf keys from the vertex descriptor of one of its 
critical vertices to navigate to an ancestor node of the key. 

More precisely, each simplex belongs to a unique quadtree leaf node and the 
simplex is conceptually assigned the key of this leaf node. In the mesh representation of the 
surface, a vertex is shared by a list of simplices (these simplices may not be instantiated). 

15 Thus, the vertex inherits a list of leaf keys from the keys assigned to the simplices that 
reference the vertex. This list, which is the vertex descriptor, is defined for all vertices in 
the surface and for the critical vertices is the same list as the one described in the paragraph 
above. 

It is not necessary to store the quadtree leaf key in a simplex because it can be 
20 computed from the vertex descriptors of the simplice's vertices, as described below. 
However, for efficiency, it can be cached at the simplex. 

Within a cell, denote by Simps(v) the set of simplices that connect to the vertex 
v. Denote by Verts(S) the set of vertices forming the comers of the simplex S. The "vertex 
descriptor" for the vertex v, is the list of leaf keys of the simplices connected to the vertex, 
25 V, . 
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Keysiy) = \J Key{Leaf{S)) 

SGSimpsiv) 

For performance, the tree is chosen such that the following condition is met: 

Card(Keysiv)) < 4 

By assigning the vertex descriptors to the vertices, it is possible to indirectly determine the 
leaf key of a simplex by 

KeyiLeaf(S)) = p| Keys(v) 

yeVertsiS) 

Having assigned the vertex descriptors to the vertices it is now a cheap operation 
to determine the criticality of a vertex at a particular depth of the tree. 
Let dim(v) = macrO'dimensipn(v), That is, let 

Q4f vis identified witha 0 — cell vertex 
dim(v)=^ 1,// visidentified witha 1 - cellvertexand noO — cellvertex 
[2,j/ visnot identified witha 1 - cell or 0 — cellvertex 

Let ki(v) be the number of ancestor keys at depth i of the vertex v, 
k. ( v) = Card ( Ancestors^ ( Keys{v)) 



10 



The vertex v is critical at level i if it satisfies any of the following: 

• The vertex v is identified with a 0-<:ell vertex, 

• The vertex v is identified with a 1-cell vertex and the number of ancestor 
keys at depth i is greater than one. 

15 • The number of ancestor keys at depth i is greater than two, 

Equivalently, the vertex v is critical at level i if, 

^.(v) >dim{v) 
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A quadtree node inherits a collection of simplices from the leaf nodes that are its 
descendants. If this collection of simplices does not reference any vertices that are 
identified with 0-cell or 1-cell vertices then, in the preferred embodiment, it is required that 
the collection of simplices must be homeomorphic to a 2-disk. In particular, the node must 
5 be connected and simply-connected (i.e. not have any holes). 

This requirement improves the efficiency of the algorithms, in particular 
migration, and saves having to perform difficult topological analysis of the quadtree nodes. 
As will be seen, this requirement will be satisfied for all algorithms of interest provided the 
initial quadtree that is built sadsfies the connectivity requirement. 
1 0 The data structures of the quadtree has the following characteristics: 

• Grid cells are assigned to quadtree leaf nodes in a regular manner, as 
described above. 

• The initial triangulation of the surface is defined by splitting grid cells. 
The simplices from a split grid-cell are then assigned to the quadtree leaf 

1 5 node to which the grid-cell was assigned. 

• Each quadtree leaf node maintains the list of simplices which it contains. 
If the simplices have not been built, this list is empty. The simplices are 
built dynamically from the grid cells when required. When the simplices 
have been built, the quadtree leaf node is marked as mesh and the grid 

2 0 representation can no longer be used. 

• Each vertex is assigned a vertex descriptor that is the list of leaf keys of 
the simplices that use the vertex. 

• Each quadtree node maintains the list of vertices that are critical for this 
quadtree node at the quadtree node's depth. 

25 • Each quadtree node that has no vertices identified to 1-cell or 0-cell ' 

vertices in any of its descendants is topologically connected and simply 
connected. 
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A quadtree node is "pure grid" if none of its descendant leaf nodes have been 

meshed. 

A "mesh" vertex is a vertex for which a simplex has been built which references 

this vertex. 

5 A "grid" vertex is a vertex for which no simplices that reference it have been 

built. 

A quadtree can be implemented in two basic forms. The first is to use a look up 
table based on the quadtree key, which can be done by using the ordering on the quadtree 
keys. The second is to maintedn links from the parent to the four children nodes similar to a 
10 linked list. The first implementation will be called a linear array quadtree, the second 
implementation a linked quadtree. 

A linear array quadtree is efficient for indexing on the key, but if large portions 
of the tree are non-existent it can be wasteful of memory. The linked quadtree is efficient 
where large portions of the tree are non-existent, but is inefficient for indexing on the key as 
15 it requires searching down the parent-child links. 

The implementation of the quadtree uses both techniques. For the coarse depths 
the linear array is used and at the finer depths a linked tree is used. This provides faster 
indexing for coarse keys, but is not too wasteful of memory. 

An iterator is provided for the tree, which maintains a key and provides efficient 
2 0 retrieval of the parent nodes. By using an iterator it is not necessary to maintain a key or the 
parent pointers in the quadtree nodes. They can be maintained in the iterator, 

A bit tree is maintained which indicates whether a particular quadtree node is 
pure grid. This allows optimization on the basis that one knows whether a particular node is 
pure grid. For example, for a surface built from a regular grid, for a pure grid node, only the 
2 5 max z and min z values need be stored to recover a bounding box. 

A bit vector for the grid cells is maintained which describes how a particular 
grid cell is split into simplices. This means it is not necessary to split all the grid cells to 
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define the simplicial structure. 

For grid vertices, there is no need to build their vertex descriptor, because keys 
can be generated on the fly for nodes that are pure grid. So vertex descriptors are only 
needed for mesh vertices. Furthermore, one knows a quadtree node that is pure grid has 
5 four neighbors and hence a pure grid node can be migrated easily. 

The implementation of the quadtree maintains the following data structures and 
implements the following algorithms. 

• A bit vector, with one bit for each grid cell, to represent the triangulation 
of grid cells. 

10 • A bit tree, with one bit for each node in the quadtree, to deiennine if a 

node is present in the tree. This is necessary to determine whether a 
node is present in persistent storage. 

• A bit tree, with one bit for each node in the quadtree, to denote whether a 
node is pure grid or not. A node can be pure grid even though its 

15 boundary vertices may connect to mesh nodes. The bit tree is an 

efficient hierarchical encoding of the hybrid grid-mesh representation. 

• An iterator is used to navigate in the tree. Thus the nodes do not need to 
store parent pointers or keys. 

• Vertex descriptors are not needed for grid vertices, 

20 • Simplices are built only when required, e.g. to support intersection 

curves. 

• The quadtree key can be implemented in 32 bits. 

• Coarse levels of the tree are implemented as a linear array quadtree. 
Finer levels are implemented as a linked quadtree. 

25 • The mesh representation stored at the quadtree leaf nodes is the micro- 

topological representation defined in the SHAPES geometry engine. 
This is to save unnecessary conversion between different mesh 
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representations when passing geometry to the SHAPES geometry 
engine. If a geometry engine had an alternate mesh representation it 
would be possible to use this representation instead. 
The implementation of the algorithms is now discussed, beginning with 
5 classification. 

At each quadtree node a bounding box is maintained which is large enough to 
contain all the simplices assigned to that quadtree node. 

To compute the intersection of two surfaces the bounding boxes can be used in a 
hierarchical manner to compute the quadtree leaf nodes that intersect. At this point, 
10 intersecting quadtree leaf nodes are resolved to their individual simplices that are then 
intersected to determine the intersection curve. 

In the preferred embodiment, the underlying geometry engine computes the 
correct topology when all possible intersecting triangles are passed to the engine. 

The next step is to make the model coherent. Making coherent consists of two 
15 fundamental steps. The first is to split simplices that lie along the macro-topological 
boundaries so the simplicial structure of the surface respects the micro-topology of its 
bounding 1 -cells and Ocelis. The second is to collect the simplices into their respective 
topologically connected components, which is called migration. 

In overview, new quadtrees are built for each new cell. As each simplex is split, 
2 0 each resulting simplex is assigned to the new quadtree for the cell. Using mesh connectivity, 
the remainder of the simplices which lie in the quadtree leaf nodes it belongs to are 
migrated. The quadtree nodes that are not split are then considered. These quadtree nodes 
are topologically connected and a flood fill using the vertex descriptors of the already 
migrated vertices is used to migrate these nodes. At this point, the critical vertices of each 
25 quadtree node in each new quadtree can be recomputed as described above. Finally, the 
bounding boxes are recomputed for each quadtree node. 

A collection of nodes, C, of the quadtree T is a "node front" if every leaf node of 
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T either has no ancestor in C or has a unique ancestor in C (a node is an ancestor of itself). 

A collection of nodes, C, of the quadtree T is a "complete node front" if C is a 
node front and every leaf node of T has an ancestor in C. 

A subtree S of T can be built from the node front C, by taking all of the 
5 ancestors of C in T. The tree structure of S is inherited from T and the collection C forms 
the leaves of S. 

The process begins with the surface being partitioned into m nodes at resolution 
level-i using an i-level subset of boundaries 288, as shown in Fig. 30a. Each level-i+1 node 
is associated with a unique level-i node 290. Each level-i node is associated with the level- 

10 i+1 nodes associated to the node 292. Each level-i node has associated with it a subset of 
the vertices that are critical at resolution level i 294. Each node at resolution level d is 
designated a leaf node 296. While the preferred embodiment has leaf nodes at the same 
level, it is also possible to have leaf nodes at numerous levels of resolution. 

Each simplex is associated with a unique leaf node 298. Each leaf node has 

15 associated with it the simplices associated to that leaf node 300. Each level-i node has 
associated with it the list of simplices which is the union of all simplices associated with the 
level-i+1 nodes grouped under the level-i node 302. The subset of boundaries for each node 
is assigned to be the boundary of the union of the simplices associated with that node 304. 
Each node is assigned a unique key and each vertex in a leaf node is assigned a key 

2 0 corresponding to that leaf node 308. 

Now assume that a second surface is classified into the model. It is determined 
which leaf nodes of the first surface intersect the leaf nodes of the second surface 310, as 
shown in Fig.30b. The intersecting simplices from the first and second surfaces are 
determined from the simplices associated to the intersecting leaf nodes 312. The original 

2 5 quadtree is split into new quadtrees and each new quadtree is associated with a new cell 
314. The subtrees of the original quadtree which have no intersecting leaf nodes are 
identified with one of the new cells 316. 
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The simplices of the first surface are split along the intersection curve 318, as 
shown in Fig. 30c. New simplices are formed by tessellating the split simplices to respect 
the macro-topology of one-cells and zero-cells passing through the original simplices 320. 
A new quadtree is built for each new cell 322. Each new simplex is assigned to the leaf 
5 node of the quadtree created for the new cell to which the new simplex belongs 324. For 
each leaf node of each new quadtree, each simplex in the original quadtree which is 
connected to a new simplex in the new quadtree leaf node and which lies in the same 
quadtree leaf node as the new simplex is migrated 326. The neighbors of a quadtree node 
are determined by finding all the keys of the critical vertices in the node 328. The coarsest 
1 0 level node which is an ancestor of a key from the critical vertices in the migrated quadtree 
nodes and which has not been split or migrated is determined and migrated to the new 
quadtree 330. 

Decimation begins with a list of critical vertices being built from the quadtree 
nodes of a complete node front 332, as shown in Fig. 30d. 

15 Figure 31 is an example of a quadtree with leaf nodes at a fixed depth. 

Examples of node fronts are {b,c,q,s,t} and {a,c,d}. Examples of complete node fronts are 
{a,b,c,d} and {a,i,j,k4,c,d}. An example of a non-node front is {a,e,b,c,d}. The collection 
{a,e,b,c,d} is not a node front because "e" does not have a unique ancestor in the collection. 

The vertices identified to one- or zero-cell vertices is removed from the list 334 

20 (Fig. 30d). All zero-cell vertices from the model which lie in the first surface are added to 
the list 336. A defined collection of one-cell vertices is added to the list 338. The 
collection of one-cell edges is recorded 340. The surface is tessellated to respect the list of 
vertices and the recorded one-cell edges 342. A requirement is imposed that the subset of 
vertices on the boundary of the first surface which are also on the boundary of the second 

2 5 surface be the same as the subset of vertices on the boundary of the second surface which 
are also on the boundary of the first surface 344. 

A first surface is partitioned into n\ nodes at resolution level-i using an i-level set " 
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of boundaries 288 (Fig. 30a). A geometrical representation of the first surface is 
maintained in persistent storage 346 (Fig. 30e). For each node, a bounding box is stored on 
a persistent storage device 348. For each node, a list of critical vertices associated with that 
node is stored on the persistent storage device 350. For each critical vertex, a vertex 
5 descriptor, a parameter value and an image value are stored on the persistent storage device 
. 352. The required portion of the first surface is loaded on demand from the persistent 
storage device 354 and that portion of the first surface not required is removed 356. A 
quadtree node is loaded on demand from persistent storage and removed when it is no 
longer needed 358. A quadtree leaf node is loaded on»demand from persistent storage and 
1 0 removed when it is no longer needed 360. 

In this way memory usage is conserved, and furthermore changes to the model 
are limited to the particular collection of sub-volumes specified. 

Even with the ability to selectively load geometry at the macro-topological level 

15 there is still a memory usage problem. For example, the user of the GQI may want to load 
the whole earth model and view it. Using the multi-resolution hierarchy it is possible to 
selectively load portions of the geometry of the surfaces. For a user who is viewing the 
whole earth model it is not necessary to load the fine details of the model. It is sufficient to 
load a collection of nodes from the coarse levels of the quadtrees to give a good 

2 0 approximation to the earth model. For a user who is viewing a small part of the whole earth 
model, it is not necessary to load the finer levels of the quadtrees not in the viewing volume. 

Consider computing the intersection of two surfaces. This requires loading all 
pairwise intersecting leaf nodes and then loading the simplices which these leaf nodes 
contain. This can be done in a recursive manner as follows. Load the bounding box for the 

2 5 root node of each quadtree of each surface. If there is no intersection between the bounding 
boxes then the two surfaces do liot intersect. If there is an intersection choose one of the 
nodes and load its children's bounding boxes and mark the node*s bounding box for 
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removal from memory. Now intersect the children's bounding boxes with the bounding box 
of the node from the other tree. If there is an intersection then recurse down the branches of 
the trees by loading the bounding boxes of thb children nodes and intersecting with the 
bounding boxes from the other tree until the leaf nodes are reached. If there is no 
5 intersection then stop the recursion. This results in pairs of intersecting quadtree leaf nodes. 
Now load all the simplices which are contained in the leaf nodes and pass them to the 
geometry engine. 

The partial loading algorithm described above for the intersection algorithm can 
be applied to all algorithms which reference infomiation in the hierarchical surface. In 
10 particular, it can be applied to the migration algorithm where nodes, simplices and vertices 
can be selectively loaded from disk. 

To achieve optimum performance, it is essential, when local changes are made 
to the model, that it is possible to map these changes to local updates in the persistent 
storage. An example of a local change is to modify the (x,y,z) position of a vertex. To 
15 maintain efficient persistent storage, it is essential this (x,y,z) position is at a limited number 
of locations in persistent storage, and preferably should be at a unique location. The 
mapping architecture is described below and is achieved by mapping the quadtree structure 
to persistent storage. 

A quadtree node consists of a fixed size component, the bounding box, and a 
2 0 variable sized component, the list of critical vertices. For efficiency, these two components 
are maintained in separate locations in persistent storage. In persistent storage the fixed size 
component models the situation in memory and uses a mix of linear array indexing with 
linked-list indexing at the finer levels. 

The architecture for storing the critical vertices is more complicated and is 
2 5 described below. 

When saving a vertex to disk it is necessary to store the vertex descriptor (the 
list of leaf keys), the parameter value of the vertex, and the image value of the vertex. 
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All vertices that are identified with 0-cell or 1-cell vertices are stored separately 
' and are loaded whenever the surface is loaded into core memory. The SHAPES geonietry 
engine maintains an index that is used to refer to the 0-cell or 1-cell vertices this vertex is 
identified with. For the geometry engine to rebuild the identification structures the vertex 
5 must store its index and on loading must restore the index. 

All other critical vertices are stored at the quadtree node, which is the finest 
common ancestor of all the leaf keys in the vertex descriptor. 

A quadtree node has three types of critical vertices, illustrated in Fig. 32. The 
type of critical vertex can be determined from the number of unique keys at the parent level. 
10 Let V be a critical vertex at level i+1. As discussed above, this implies: 

*w(v)>dim(v) 

Accordingly, v must be either parent critical, edge critical or sub-critical: 

1. V is "parent critical"(e.g. node 362, Fig. 32) if: 

/:,.(v)>dim(v) 

2. V is "edge critical" (e.g. nodes 364 and 366) if the vertex was not critical at 
the parent level and lies on the boundary of the parent node: 

A:, (v) = dim(v) 

15 3. V is "sub-critical" (e,g. node 368) if the vertex is neither parent critical nor 

edge critical: 

/:,.(v)<dim(v) 

The critical vertices of a quadtree node are maintained in three lists, the parent 
critical, the edge critical and the sub-critical vertices. For the parent critical, vertices the list 
is a list of indices into the parent critical vertices. For the sub-critical vertices, the vertex is 
2 0 stored at the parent node, because the parent is the finest common ancestor of the list of leaf 
keys. The list of sub-critical vertices is a list of indices into the vertices stored at the parent 
node. Identifying an edge critical vertex requires a quadtree key, which specifies the node, 
which stores the vertex, together with an index into the list of stored vertices of the node 
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that stores the vertex. 

The list of vertices which is stored at a particular node can be broken into a 
collection of buckets. For example, every vertex has a unique depth, this being the depth 
when the vertex first becomes critical. For the quadtree nodes that arc at the vertex's depth 
5 and have the vertex as a critical vertex, the vertex must be either edge critical or sub-critical 
for all these quadtree nodes. Thus, the list of stored vertices can be broken into vertices that 
will be edge critical vertices and vertices that will be sub-critical vertices. The edge critical 
vertices can be broken down further by finding the pair of quadtree keys at the depth of the 
vertex and by identifying the child index of each key in its parent which gives a pair of two 

10 bit indices which when ordered can be used to identify a bucket. The edge bucket can be 
fiirther broken down by using the depth of each vertex. In fact, an edge bucket inherits a 
binary tree ftom the quadtree structure (it is the restriction of the quadtree to the particular 
edge), and this can be used to provide hierarchical loading of the edge critical vertices. 

The quadtree leaf node has an additional component beyond the other nodes in 

15 the tree and that is the mesh representation of the node. This is a simple list of triples of 
indices specifying the indices into the critical vertices, edge vertices and internal vertices of 
this quadtree node. As illustrated in Fig. 33, a grid representation of a surface is stored, the 
grid being made up of grid cells 370. A mesh representation of a portion of the surface is 
then formed by triangulating a subset of the grid cells 372. 

20 Since the quadtree leaf node must contain the list of simplices, the mesh content 

of the leaves is stored in a separate location on the disk. This storage is implemented using 
a blocked linked list. This allows maximum flexibility for adding and deleting simplices 
from the database, but provides an efficient means to load a particular leaf node. 

It is also possible to implement a multiresolution hierarchy for multivalued 

2 5 surfaces, as discussed below. 

The GQI must be able to import a triangulated surface defined as an Inventor 
face set. The simplices in the face set can be oriented in a completely arbitrary manner, and * 
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overall the surface can be multi-valued and non-manifold. While it is possible to formally 
apply the quadtree partitioning scheme to such a surface, it is unlikely that the tiling element 
patches will be connected. What is needed is a mapping of the surface onto a planar region 
that preserves the surface's simplicial connectivity. 
5 One approach for constructing this mapping is to model the surface deformation 

as an energy minimization problem whose solution is consistent with Hamilton's Least 
Action Principle: 

Given an elastic surface g(u,v,t) = (x(u,v,t), y(u.v,t), z(u,v,t)) that is to be 
deformed into a shape g'(u,v,t) = (x'(u,v,t), y'(u,v,t), z'(u,v,t)). Assume that a set 
10 of imaginary massless springs join sample points on g to their final place on g'. 

Then the motion of deformation follows a path in time such that the action 

b 
a 

is minimized. (The time interval is [a,b], K is the surface system's kinetic 
energy, and P is its potential energy. The difference K-P is the system's 
"LaGrangian".) 

15 In order to apply this principle to triangulated surfaces, the surface must be "sampled" at its 
vertices. An application can constrain the potential energy in time three ways: 

1. By the position of each vertex. The position constraint s^plied to the 
potential energy stored in the springs is defined by 

\ljC(ig'-g)'dudv 
where C is a spring constant that is a function of the (u,v) sample 
20 coordinate. Frequently, C is independent of the (u,v) sample coordinate. 

The Euler-LaGrange theorem applied to the minimization of this 
integral says that the deformation gXu,v,t) must converge to g(u,v,t) over 
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time. 



2. By the orientation of the surface at each sample vertex. The orientation 
constraint is defined by 



where D denotes a different spring constant. Again, D is usually set to a 
fixed common value. The Euler-LaGrange theorem applied to the minimizing 
function for this integral says that the deformation g(u,v,t) must be a harmonic 
function. This is good, because a harmonic function causes the least distortion 
to the aspect ratio of the simplices on the deformed surface. A sliver results 
when the aspect ratio of a simplex is poor. Therefore, if the intention is to 
tessellate the planar deformation, then minimizing the number of slivers and 
their degree of malformation is important. 

3. Finally, the curvature constraint applied to the spring's potential energy 
is given by a minimizing function for the integral equation. 



where E is a spring constant. Again, E is usually set to a fixed coimnon 
value. The Euler-LaGrange theorem says that a minimizing function for 
this constraint satisfies a biharmonic relation. This form of constraint is 
used in some grid-based surface modeling systems. Typically, the fiill 
svirface is subdivided into small patches. On each patch a solution is 
computed, and some form of smoothing is applied across the interface 
between two patches. 
A GQI surface is triangulated, so the minimizing function must preserve the 
of every vertex. Hence a linear harmonic deformation is ideal. Hamilton's 
e says nothing about how vertex pairs, i.e., the edges that form the topology,- 
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deform. This is significant, because a smooth vertex deformation can reposition the vertex 
endpoints of a set of connected simplex edges so that the deformed edges cross, failing to 
preserve the surface's topological definition. When this happens, the GQI has to split the 
surface before the folding occurs and restart the process on the remainder of the surface. 
5 Thus, the deformation of a surface feature may be defined as a large number of cells. A 
large cell count slows down classification, and forces the GQI to create a quadtree for each 
2-cell, exploding memory. What is needed is a method to parameterize the surface given 
just the surface's vertex and simplex connectivity, i.e., a topology based method. 

A method is described below for constmction of a single 2-cell for the case that 

10 the surface is oriented with an arbitrary number of holes. This method does not control 
metric distortion, in contrast to a harmonic deformation. The GQI computes this 
deformation in order to construct a quadtree rather than as preparation for further 
tessellation, so some distortion is acceptable. When applied to a gridded triangulated 
surface, the mapping is equivalent to the projection of the surface onto one of the coordinate 

15 planes. 

Let F, E, and V be the number of faces, edges, and vertices, respectively, of a 
triangulated surface. The Euler Characteristic X of the surface is defined asX = F- E + V. 
When the surface is presented in "normal" form, the Euler Characteristic X of an oriented 
surface is equal to X = 2 - (2 • g + r), where g denotes the genus of the surface (which 
2 0 equals the number of embedded tori) and r is the number of boundary curves ("holes"), 
grouped into GQI frames. Fig. 34 lists some common surfaces and their Euler 
Characteristics. 

A standard result in topology is that two oriented surfaces are homeomorphic if 
and only if they have the same genus and number of frames. The Euler Characteristic is a 
2 5 function of genus and the number of frames, so it follows that the Euler Characteristic is 
independent of the triangulation. The GQI supports queries for the evaluation of F, E, and 
V, as well as frames. Thus two oriented triangulated surfaces that have the same Euler 



wo 00/19380 PCTAJS99/22194 



-42- 

Characteristic and the same number of frames must be homeomorphic. 

In geological applications, multivalued surfaces such as a recumbent fault, a salt 
body, or a horizontal wellbore are encountered. All three surfaces have genus zero and have 
Euler Characteristic 0, 1 or 2, assuming that the surface has zero internal holes. The 
5 quadtree construction assumes that the surface is single- valued in its (u,v) coordinates, so to 
apply the quadtree construction to a multi-valued surface, it is necessary to construct a 
homeomorphism between the surface and a single-valued representation with the correct 
Euler Characteristic. If the original surface has Euler Characteristic 1, e.g., a recumbent 
fault, a good choice is to map it onto a planar square or family of connected squares. If the 

10 original surface has Euler Characteristic 0, e.g., a wellbore, then the square should be 
replaced by a squared annulus. Finally if the surface has Euler Characteristic 2, then it can 
be split into two equal-area parts with each part of Euler Characteristic 1 with the previous 
construction applied to each part. 

Here are the steps in the construction of the homeomorphism between a set of 

15 connected planar squares and an oriented genus 0 triangulated surface with Euler 
Characteristic 

1. Build the surface as a non-parametric web, so that the non-manifold 
boundaries are apparent. From now on, assume that the surface is a 2- 
manifold. 

20 2. Find the surface's external boundary curve, which is defined by the 

zeroth GQI frame instance. Analyze the external boundary for 
embedded disc obstructions, as shown in Fig. 35. Remove each Type (a) 
obstruction by adding a new vertex plus edges to form a new simplex 
that uses the obstructing vertex as a comer, as shown in Fig. 36. If a 

25 Type (b) obstmction is detected, then exclude it from the remainder of 

the algorithm and at the conclusion add it to the final planar square. If a 
Type (c) obstruction is detected, then split the surface into three parts. 
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building two squares and a connector, essentially reducing to a Type (b) 
problem. 

Compute the area of the (possibly extended) surface, ignoring Type (b) 
obstructions, and construct a square of that size. The area of each surface 
simplex can be computed from the cross product of two sides of the 
simplex, thought of as vectors. Randomly select four points along the 
external boundary that are spaced roughly 1/4 of the contour apart from 
each other. These points are mapped to the comers of the square. The 
remaining points are assigned in relative proportion to the appropriate 
side of the square. Initialization is now complete. 
Determine the contour formed from the set of all vertices connected to 
the most recently defined contour, as shown in Fig. 37. First find four 
vertices on the new contour that are connected to the comers of the old 
contour. If no edge exists between a comer to the new contour, then 
split the obstructing simplex. Note that this might cause other simplices 
to split. A good choice of a second vertex is the simplex 's barycenter. In 
part (a) of Fig. 37, the dotted edge has been added to join the two 
comers. 

If two edges of the contour form a Type (a) obstruction, then redefine the 
new contour to include that simplex unless the obstruction simplex's 
valence 2 vertex is directly connected to a comer vertex of the previous 
contour. In part (b) of Fig. 37, the dark simplex is a Type (a) 
obstruction. Assuming that the valence 2 vertex is not directly 
connected to a comer vertex of the previous contour, remove the 
obstruction by incorporating the obstructing simplex 's remaining edge 
into the new contour, effectively flattening the contour. In Fig. 37 the 
two offending interior comers do not obstruct the direct connection to " 
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comer vertices of the previous contour, so both offenders can be safely 
ignored. 

6. Since the full surface is oriented and of genus 0, the set of simplices 
connecting the two contours is also oriented and of genus 0. Assuming 

5 that the surface has no holes and contains no Type (a) obstmctions and 

no non-manifold linkage of the interior and exterior contours, it follows 
that the surface patch defined by this set of simplices has Euler 
Characteristic. Referring to Fig. 38, the patch is homeomorphic to the 
lateral surface of a tmncated square-based pyramid. Deform the surface 

1 0 formed by the simplices between the two contours into the lateral surface 

of a truncated square-based pjrramid, then project the upper contour onto 
the interior of the square defined by the lower contour, as shown in Fig. 
38. Note that two edges on a lateral face cross if and only if their 
projection cross, explaining why the choice of a lateral surface of a 

15 square-based pyramid. If a comer simplex was split in step #4, then 

replace its split parts by the projection of the original simplex. 

7. It is possible that the interior and exterior contours together describe a 
connected sum of lateral surfaces of square-based pyramids, due to Type 
(a) obstmctions, or meet in a non-manifold manner, see Fig. 39. An 

2 0 example of this would be an otherwise flat surface on which two hills 

exist. If this happens, the process divides the surface into the number of 
connected components identified by the interior contour and applies the 
construction to each component. (Note that the hexagonal region 
between the two shaded areas has Euler Characteristic 1, so it is 

25 homeomorphic to a disc.) When this happens, the algorithm is applied 

to each sub-square and any Type (a) obstmctions are fiJled in. A 
homeomorphism preserves continuity, so the number of embedded 
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squares equals the number of hills and valleys in the multi-valued 
surface. Note that running the construction in reverse is a simple way to 
create a surface with any number of hills and valleys, salt domes, etc. 

8. Repeat this identification until the current contour is the boundary of an 
5 embedded disk, as shown in Fig. 40. Complete the square construction 

by deforming the indicated disk to a square and adding it in. 

9. It is possible that there may not be four distinct vertices on a contour, 
e.g., a pyramid formed with a triangle cap and a square base. If this 
happens, then the nested square is replaced by a triangle and the 

i 0 construction proceeds as above, as shown in Fig. 41. 

Loosely speaking, the homeomorphism "melts" the surface onto a plane. 
Consequently, points close together in (x,y,z) space may not be close together in the plane. 
(This is another way of saying that this method does not control metric distortion.) An 
extreme example is a spherical surface with a tiny hole at a pole. If the hole is used as the 

15 initial boundary, then the triangles forming the boundary can be quite distorted in the 
mapping plane. When the surface is single-valued, it may be acceptable to use the 
projection of the bounding box onta the appropriate plane as the initial square and create 
nested square-sided armuli that proportional with respect to area. In any case, distortion is 
not a problem, since the planar version of the surface is used only to facilitate the quadtree 

20 decomposition of the surface. The GQI constructs the quadtree on the homeomorphic 
planar region, then lifts it back. The only significant change in the quadtree hierarchy API 
is that ray picking can return more than one tiling element even when the ray does not pick a 
sirhplex comer. 

This construction can be applied to a surface with an arbitrary number of holes 
25 in it by adding to the area of the tiling cover of each hole to the area of the original surface. 
Also, this approach can tile a toroidal shape. Indeed, all that is needed is to randomly 
choose two closed paths in opposing directions, then "unroll" the toroidal shape along the 
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two paths. Assuming that the surface has no holes, the two formulas for the Euler 
Characteristic can be equated to infer when a toroidal shape is present and the quadtree 
construction applied to each toroidal element. Not much is made of this, however, because 
toroidal structures are seldom encountered in geological modeling. 
5 The invention has application outside the field of geological modeling. In 

particular, the invention has application in any field in which data are presented 
geometrically and graphically at more than one level of resolution. For example, the 
invention would be useful in representing medical models, such as those developed in 
magnetic resonance imaging ("MRT'). Further, the invention would be useful in 

1 0 representing computer aided design ("CAD") models. 

The invention may be implemented in hardware or software, or a combination of 
both. However, preferably, the invention is implemented in computer programs executing 
on programmable computers each comprising a processor, a data storage system (including 
volatOe and non-volatile memory and/or storage elements), at least one input device, and at 

15 least one output device. Program code is applied to input data to perform the functions 
described above and generate output information. The output information is applied to one 
or more output devices, in known fashion. 

Each program is preferably implemented in a high level procedural or object 
oriented programming language (such as C++ or C) to communicate with a computer 

2 0 system. However, the programs can be implemented in assembly or machine language, if 
desired. In any case, the language may be a compiled or an interpreted language. 

Each such computer program is preferably stored on a storage media or device 
(e.g., ROM or magnetic/optical disk or diskette) readable by a general or special purpose 
programmable computer, for configuring and operating the computer when the storage 

2 5 media or device is read by the computer to perform the procedures described herein. The 
inventive system may also be considered to be implemented as a computer-readable storage 
medium, configured with a computer program, where the storage medium so configured 
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causes a computer to operate in a specific and predefined manner to perform the functions 
described herein. 

Other embodiments are within the scope of the following claims. 

What is claimed is: 
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1 LA method for representing a first surface at multiple levels of 

2 resolution, the first surface comprising zero or more zero-cells, zero or more one- 

3 cells and one or more two cells, the method being implemented in a programmed 

4 computer comprising a processor, a memory, a persistent storage system, at least one 

5 input device, and at least one output device, the method and a model being stored on 

6 a computer-readable media, the method representing the model on one of the output 

7 devices, the method comprising: 

8 partitioning the first surface with one or more boundaries, each level of 

9 resolution having a subset of the boundaries. 

1 2. The method of claim 1 fiirther comprising: 

2 partitioning the first surface into ni nodes at resolution level-i using the level- 

3 i subset of boundaries; 

4 associating each level-i+1 node with a unique level-i node; 

5 associating with each level-i node the level-i+1 nodes associated to the node. 

1 3. The method of claim 2 ftirther comprising: 

2 associating with each level-i node a subset of vertices that are critical at 

3 resolution level i. 

1 4. The method of claim 3, wherein level d is the deepest level of 

2 resolution and wherein the first surface is divided into simplices. further comprising: 

3 designating each node at resolution level d a leaf node; 

4 associating each simplex to a unique leaf node; and 

5 associating with each leaf node the simplices associated to that leaf node. 
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1 5. The method of claim 4 further comprising: 

2 associating with a level-i node the list of simplices which is the union of all 

3 simplices associated with the level-i+1 nodes associated to the level-i 

4 node. 

1 6. The method of claim 5 further comprising assigning the subset of 

2 boundaries for each node the boundary of the union of the simplices associated with 

3 that node. 

1 7. The method of claim 6 wherein the nodes form an original tree and each 

2 node is assigned a unique key. 

1 . 8. The method of claim 7 further comprising: 

2 assigning to each vertex in a leaf node the key corresponding to that leaf 

3 node. 

1 9. The method of claim 8 further comprising: 

2 storing the representation of a second surface in the computer-readable 

3 media, the second surface having nodes, leaf nodes, vertices, critical 

4 vertices and simplices. 

1 10. The method of claim 9 further comprising: 

2 determining which leaf nodes of the first surface intersect the leaf nodes of 

3 the second surface; 

4 determining the intersecting simplices from the first and second surfaces 

5 from the simplices associated to the intersecting leaf nodes. 

1 11. The method of claim 4 wherein each node except the leaf nodes has a 

2 subtree, further comprising splitting the original tree into new trees and associating 

3 each new tree with a new cell. 
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12. The method of claim 11 further comprising identifying the subtrees of 
the original tree which have no intersecting leaf nodes with one of the new cells. 

13. The method of claim 12 further comprising: 

splitting the simplices of the first surface along the intersection curve; 
forming new simplices by tessellating the split simplices to respect the 

macro-topology of one-cells and zero-cells passing through the 

original simplices; 
building a new tree for each new cell; 

assigning each new simplex to the leaf node of the tree created for the new 
cell to the which the new simplex belongs; 

for each leaf node of each new tree, migrating each simplex in the original 
tree which is connected to a new simplex in the new tree leaf node 
and which lies in the same tree leaf node as the new simplex; 

determining the neighbors of a tree node by fmding all the keys of all the 
critical vertices in the node; and 

determining the coarsest level node which is an ancestor of a key from the 
critical vertices in the migrated tree nodes and has not been split or 
migrated and migrating that node to the new tree. 



& 
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1 14. The method of claim 8, wherein a complete node front of the tree of the 

2 first surface and a collection of vertices on a boundary of the first surface are defined, 

3 further comprising: 

4 building a list of critical vertices from the tree nodes of the complete node 

5 front; 

6 removing from the list those vertices identified to one- or zero-cell vertices; 

7 adding to the list all zero-cell vertices from the model which lie in the first 

8 surface; 

9 adding to the list the defined collection of one-cell vertices; 

10 recording the collection of one-cell edges; and 

11 tessellating the surface to respect the list of vertices and the recorded one-cell 

12 edges. 

1 15. The method of claim 14 further comprising: 

2 requiring the subset of vertices on the boundary of the first surface which are 

3 also on the boundary of the second surface to be the same as the 

4 subset of vertices on the boundary of the second surface which are 

5 also on the boundary of the first surface. 

1 16. The method of claim 1 further comprising: 

2 maintaining in a persistent storage a geometrical representation of the first 

3 surface. 

1 17. The method of claim 2 further comprising: 

2 storing for each node a bounding box on a persistent storage device, and 

3 storing for each node a list of critical vertices associated with that 

4 node on the persistent storage device. 
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18. The method of claim 17 wherein storing the list of critical vertices 
comprises: 

storing a vertex descriptor for each vertex; 
storing a parameter value for each vertex; and 
storing an image value for each vertex. 

19. The method of claim 18 further comprising loading on demand from 
persistent storage into memory that portion of the first surface required and removing 
from memory that portion of the first surface not required. 

20. The method of claim 19 wherein a tree node is loaded from persistent 
storage on demand and is removed from memory when no longer required. 

21. The method of claim 20 wherein the simplices associated with a tree leaf 
node are loaded from persistent storage on demand and removed from memory when 
no longer required. 

22. A method for representing a first surface at multiple levels of resolution, 
the method being implemented in a programmed computer comprising a processor, a 
data storage system, at least one input device, and at least one output device, the 
method and a model being stored on a computer-readable media, the method 
representing the model on one of the output devices, the method comprising: 

storing a grid representation of the first surface, the grid representation being 

made up of grid cells; 
forming a mesh representation of a portion of the first surface by 

triangulating a subset of the grid cells; and 
inserting the first surface into the model. 
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1 23. A computer system for representing a first surface at multiple levels of 

2 resolution, the first surface comprising zero or more zero-cells, zero or more one- 

3 cells and one or more two cells, the computer system comprising a processor, a data 

4 storage system, at least one input device, and at least one output device, the first 

5 surface being stored on the data storage system, the computer system comprising: 

6 means for partitioning the first surface with one or more boundaries, each 

7 level of resolution having a subset of the boundaries. 

1 24. The computer system of claim 23 further comprising: 

2 means for partitioning the first surface into n\ nodes at resolution level-i 

3 using the level-i subset of boundaries; 

4 means for associating each level-i+1 node with a unique level-i node; 

5 means for associating with each level-i node the level-i+1 nodes associated to 

6 the node. 

1 25. The computer system of claim 24 further comprising: 

2 means for associating with each level-i node a subset of vertices that are 

3 critical at resolution level i. 

1 26. The computer system of claim 25, wherein level d is the deepest level of 

2 resolution and wherein the first surface is divided into simplices, further comprising: 

3 means for designating each node at resolution level d a leaf node; 

4 means for associating each simplex to a unique leaf node; and 

5 . means for associating with each leaf node the simplices associated to that leaf 

6 node. 

1 27. The computer system of claim 26 further comprising: 

2 means for associating with a level-i node the list of sirnplices which is the 

3 union of all simplices associated with the level-i+1 nodes associated 

4 to the level-i node. 
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1 28. The computer system of claim 27 further comprising means for assigning 

2 the subset of boundaries for each node the boundary of the union of the simplices 

3 associated with that node. 

1 29. The computer system of claim 28 wherein the nodes form an original tree 

2 and each node is assigned a unique key. 

1 30. The computer system of claim 29 further comprising: 

2 means for assigning to each vertex in a leaf node the key corresponding to 

3 that leaf node. 

1 31. The computer system of claim 30 further comprising: 

2 means for storing the representation of a second surface in the computer- 

3 readable media, the second surface having nodes, leaf nodes, vertices, 

4 critical vertices and simplices. 

1 32. The computer system of claim 31 further comprising: 

2 means for determining which leaf nodes of the first surface intersect the leaf 

3 nodes of the second surface; 

4 means for determining the intersecting simplices from the first and second 

5 surfaces from the simplices associated to the intersecting leaf nodes. 

1 . 33. The computer system of claim 26 wherein each node except the leaf 

2 nodes has a subtree, further comprising means for splitting the original tree into new 

3 trees and associating each new tree with a new cell. 

1 34. The computer system of claim 33 further comprising means for 

2 identifying the subtrees of the original tree which have no intersecting leaf nodes 

3 with one of the new cells. 
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1 35. The computer system of claim 34 further comprising: 

2 means for splitting the simplices of the first surface along the intersection 

3 curve; 

4 means for forming new simplices by tessellating the split simplices to respect 

5 the macro-topology of one-cells and zero-cells passing through the 

6 original simplices; 

7 means for building a new tree for each new cell; 

8 means for assigning each new simplex to the leaf node of the tree created for 

9 the new cell to the which the new simplex belongs; 

10 for each leaf node of each new tree, means for migrating each simplex in the 

11 original tree which is connected to a new simplex in the new tree leaf 

12 node and which lies in the same tree leaf node as the new simplex; 

13 means for determining the neighbors of a tree node by finding all the keys of 

14 all the critical vertices in the node; and 

15 means for determining the coarsest level node which is an ancestor of a key 

1 6 from the critical vertices in the migrated tree nodes and has not been 

1 7 split or migrated and migrating that node to the new tree. 
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36. The computer system of claim 30. wherein a complete node front of the 
\ of the first surface and a collection of vertices on a boundary of the first surface 
defined, further comprising: 

means for building a list of critical vertices from the tree nodes of the 

complete node front; 
means for removing from the list those vertices identified to one- or zero-cell 

vertices; 

means for adding to the list all zero-cell vertices from the model which lie in 
the first surface; 

means for adding to the list the defined collection of one-cell vertices; 
means for recording the collection of one-cell edges; and 
means for tessellating the surface to respect the list of vertices and the 
recorded one-cell edges. 

37. The computer system of claim 36 further comprising: 

means for requiring the subset of vertices on the boundary of the first surface 
which are also on the boundary of the second surface to be the same 
as the subset of vertices on the boundary of the second surface which 
are also on the boundary of the first surface. 

38. The computer system of claim 23 further comprising: 

means for maintaining in a persistent storage a geometrical representation of 
the first surface. 

39. The computer system of claim 24 further comprising: 

means for storing for each node a bounding box on a persistent storage 
device, and storing for each node a list of critical vertices associated 
with that node on the persistent storage device. 
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1 40. The computer system of claim 39 wherein the means for storing the list 

2 of critical vertices comprises: 

3 means for storing a vertex descriptor for each vertex, 

4 means for storing a parameter value for each vertex, and 

5 means for storing an image value for each vertex. 

1 41. The computer system of claim 40 further comprising means for loading 

2 on demand from persistent storage into memory that portion of the first surface 

3 required and means for removing from memory that portion of the first surface not 

4 required. 

1 42. The computer system of claim 41 wherein the means for loading loads a 

2 tree node from persistent storage on demand and the means for removing removes a 

3 tree node from memory when it is no longer required. 

1 43, The computer system of claim 42 wherein the means for loading loads 

2 simplices associated with a tree leaf node from persistent storage on demand and the 

3 means for removing removes simplices from memory when they are no longer 

4 required. - 

1 44. A computer system for representing a first surface at multiple levels of 

2 resolution, the computer system comprising a processor, a data storage system, at 

3 least one input device, and at least one output device, the first surface being stored on 

4 the data storage system, the computer system comprising: 

5 means for storing a grid representation of the first surface, the grid 

6 representation being made up of grid cells; 

7 means for forming a mesh representation of a portion of the first surface by 

8 triangulating a subset of the grid cells; and 

9 means for inserting the first surface into the model. 
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1 45. A computer program, residing on a computer-readable medium, 

2 comprising instructions for causing a computer, comprising a processor, a data 

3 storage system, at least one input device, and at least one output device, to represent 

4 a first surface at multiple levels of resolution, the first surface comprising zero or 

5 more zero-cells, zero or more one-cells and one or more two cells, the computer 

6 program comprising instructions for causing the computer to: 

7 partition the first surface with one or more boundaries, each level of 

8 resolution having a subset of the boundaries. 

1 46, The computer program of claim 45 further comprising instructions for 

2 causing the computer to: 

3 partition the first surface into m nodes at resolution level-i using the level-i 

4 subset of boundaries ; 

5 associate each level-i+1 node with a unique level-i node; 

6 associate with each level-i node the level-i+1 nodes associated to the node. 

1 47. The computer program of claim 46 further comprising instructions for 

2 causing the computer to: 

3 associate with each level-i node a subset of vertices that are critical at 

4 resolution level i. 

1 48. The computer program of claim 47, wherein level d is the deepest level 

2 of resolution and wherein the first surface is divided into simplices, further 

3 comprising instructions for causing the computer to: 

4 designate each node at resolution level d a leaf node; 

5 associate each simplex to a unique leaf node; and 

6 associate with each leaf node the simplices associated to that leaf node. 
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1 49. The computer program of claim 48 further comprising instnactions for 

2 causing the computer to: 

3 associate with a ievel-i node the list of simplices which is the union of all 

4 simplices associated with the level-i+1 nodes associated to the level-i 

5 node. 

1 50. The computer program of claim 49 further comprising instructions for 

2 causing the computer to assign the subset of boundaries for each node the boundary 

3 of the union of the simplices associated with that node. 

1 51. The computer program of claim 50 wherein the nodes form an original 

2 tree and each node is assigned a unique key. 

1 52. The computer program of claim 51 further comprising instructions for 

2 causing the computer to: 

3 assign to each vertex in a leaf node the key corresponding to that leaf node. 

1 53. The computer program of claim 52 further comprising instructions for 

2 causing the computer to: 

3 store the representation of a second surface in the computer-readable media, 

4 the second surface having nodes, leaf nodes, vertices, critical vertices 

5 and simplices. 

1 54. The computer program of claim 53 further comprising instructions for 

2 causing the computer to: 

3 determine which leaf nodes of the first surface intersect the leaf nodes of the 

4 second surface; and 

5 determine the intersecting simplices from the first and second surfaces from 

6 the simplices associated to the intersecting leaf nodes. 
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1 55. The computer program of claim 48 wherein each node except the leaf 

2 nodes has a subtree, further comprising instructions for causing the computer to split 

3 the original tree into new trees and to associate each new tree with a new cell. 

1 56. The computer program of claim 55 further comprising instructions for 

2 causing the computer to identify the subtrees of the original tree which have no 

3 intersecting leaf nodes with one of the new cells. 

1 57. The computer program of claim 56 further comprising instructions for 

2 causing the computer to: 

3 split the simplices of the first surface along the intersection curve; 

4 form new simplices by tessellating the split simplices to respect the macro- 

5 topology of one-cells and zero-cells passing through the original 

6 simplices; 

7 build a new tree for each new cell; 

8 assign each new simplex to the leaf node of the tree created for the new cell 

9 to the which the new simplex belongs; 

10 for each leaf node of each new tree, migrate each simplex in the original tree 
which is connected to a new simplex in the new tree leaf node and 

3^2 which lies in the same tree leaf node as the new simplex; 

13 determine the neighbors of a tree node by finding all the keys of all the 

1 4 critical vertices in the node; and 

15 determine the coarsest level node which is an ancestor of a key from the 

16 critical vertices in the migrated tree nodes and has not been split or 

17 migrated and migrate that node to the new tree. 
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1 58. The computer program of claim 52, wherein a complete node front of the 

2 tree of the first surface and a collection of vertices on a boundary of the first surface 

3 are defined, further comprising instructions for causing the computer to: 

4 build a list of critical vertices from the tree nodes of the complete node front; 

5 remove from the list those vertices identified to one- or zero-cell vertices; 

6 add to the list all zero-cell vertices from the model which lie in the first 

7 surface; 

8 add to the list the defined collection of one-cell vertices; 

9 record the collection of one-cell edges; and 

10 tessellate the surface to respect the list of vertices and the recorded one-cell 

11 edges. 

1 59. The computer program of claim 58 further comprising instructions for 

2 causing the computer to: 

3 require the subset of vertices on the boundary of the first surface which are 

4 also on the boundary of the second surface to be the same as the 

5 subset of vertices on the boundary of the second surface which are 

6 also on the boundary of the first surface. 

1 60. The computer program of claim 45 further comprising instructions for 

2 causing the computer to: 

3 maintain in a persistent storage a geometrical representation of the first 

4 surface. 

1 61. The computer program of claim 46 further comprising instructions for 

2 causing the computer to: 

3 store for each node a bounding box on a persistent storage device, and store 

4 for each node a list of critical vertices associated with that node on the 

5 persistent storage device. 
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62. The computer program of claim 61 wherein the instructions for causing 
the computer to store the list of critical vertices comprises instructions for causing 
the computer to: 

store a vertex descriptor for each vertex, 
store a parameter value for each vertex, and 
store an image value for each vertex. 

63. The computer program of claim 62 further comprising instructions for 
causing the computer to load on demand from persistent storage into memory that 
portion of the first surface required and remove from memory that portion of the first 
surface not required. 

64. The computer program of claim 63 wherein the computer program 
causes a tree node to be loaded from persistent storage on demand and to be removed 
from memory when no longer required. 

65. The computer program of claim 64 wherein the computer program 
causes simplices associated with a tree leaf node to be loaded from persistent storage 
on demand and removed from memory when no longer required. 

66. A computer program, residing on a computer-readable medium, 
comprising instructions for causing a processor, a data storage system, at least one 
input device, and at least one output device, to represent a first surface at multiple 
levels of resolution, the computer program comprising instructions for causing the 
computer to: 

store a grid representation of the first surface, the grid representation being 

made up of grid cells; 
form a mesh representation of a portion of the first surface by triangulating a 

subset of the grid cells; and 
insert the first surface into the model. 
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1 67. The method of claim 1 further comprising parameterizing the first 

2 surface, wherein the first surface has one or more terminating contours and is 

3 tessellated into triangles, wherein parameterizing comprises 

4 if the first surface is not of Euler Characteristic 1 : 

5 cutting the first surface into patches, each patch being of Euler 

6 Characteristic 1. 
7 

1 68. The method of claim 67 wherein parameterizing comprises 

2 repeating until reaching, all terminating contours and all triangles 

3 have been processed; 

4 if an annulus of triangles around the boundary of the 

5 remaining triangles in a patch cannot be identified 

6 because of an obstruction: 

7 refining the obstruction by adding exterior edges until 

8 the obstruction is absorbed into the annulus; 

9 constructing a parameterization of the annulus. 

1 69. The method of claim 1, further comprising 

2 filling a hole in the first surface with a tile. 

1 70. The method of claim 69 further comprising parameteri2dng the first 

2 surface, wherein the first surface has one or more terminating contours and is 

3 tessellated into triangles, wherein parameterizing comprises 
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4 repeating until reaching all terminating contours and all triangles 

5 have been processed: 

6 if an annulus of triangles around the boundary of the 

7 remaining triangles of the first surface with the hole 

8 filled cannot be identified because of an obstruction: 

9 refining the obstruction by adding exterior edges until 

10 the obstruction is absorbed into the annulus; 

11 constructing a parameterization of the annulus. 

1 71. The computer system of claim 23 further comprising means for 

2 parameterizing the first surface, wherein the fu-st surface has one or more terminating 

3 contours and is tessellated into triangles, wherein the means for parameterizing 

4 comprises 

5 if the first surface is not of Euler Characteristic 1 : 

6 means for cutting the first surface into patches, each patch 

7 being of Euler Characteristic 1 . 



1 
2 



72. 

comprises 



The computer system of claim 71, wherein the means for parameterizing 
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3 means for repeating until reaching all teiminating contours and all 

4 triangles have been processed: 

5 if an annulus of triangles around the boundary of the 

6 remaining triangles in a patch cannot be identified 

7 because of an obstruction: 

8 means for refining the obstruction by adding exterior 

9 edges until the obstruction is absorbed into the 

10 annulus; 

11 means for constructing a parameterization of the annulus. 

1 73. Thecomputer system of claim 23, further comprising 

2 means for filling a hole in the first surface with a tile. 

1 74. The computer system of claim 73, wherein the means for parameterizing 

2 comprises 

3 means for repeating until reaching all terminating contours and all 

4 triangles have been processed: 

5 if an annulus of triangles around the boundary of the 

6 remaining triangles of the first surface with the hole 

7 filled cannot be identified because of an obstruction: 

8 means for refining the obstruction by adding exterior 

9 edges until the obstruction is absorbed into the 

10 annulus; 

1 1 means for constructing a parameterization of the 

12 annulus. 
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75. The computer program of claim 45 further comprising instructions for 
causing the computer to parameterize the first surface, wherein the first surface has 
one or more terminating contours and is tessellated into triangles, wherein the 
computer program further comprises 

if the first surface is not of Euler Characteristic 1, instructions for 
causing the computer to: 

cut the first surface into patches, each patch being of Euler 
Characteristic 1. 

76. The computer program of claim 75 wherein the instructions for causing 
the computer to parameterize the first surface comprise 

instructions for causing the computer to repeat, until reaching all 
terminating contours and all triangles have been processed: 
if an annulus of triangles around the boundary of the 

remaining triangles in a patch cannot be identified 
because of an obstruction: 

instructions for causing the computer to refine the 

obstruction by adding exterior edges until the 
obstruction is absorbed into the annulus; 
instructions for causing the computer to construct a 
parameterization of the annulus. 

77. The computer program of claim 45, further comprising instructions for 

causing the computer to 

fill a hole in the first surface with a tile. 
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1 78. The computer program of claim 77 wherein the instructions for causing the 

2 computer to parameterize the first surface comprise 

3 instructions for causing the computer to repeat, until reaching ail 

4 terminating contours and all triangles have been processed: 

5 if an annulus of triangles around the boundary of the 

6 remaining triangles of the first surface with the hole 

7 filled cannot be identified because of an obstruction: 

8 instructions for causing the computer to refine the 

9 obstruction by adding exterior edges until the 

1 0 obstruction is absorbed into the annulus; 

11 instructions for causing the computer to construct a 

12 parameterization of the annulus. 
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